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We derive the most general junction conditions for the fourth-order brane gravity constructed of 
arbitrary functions of curvature invariants. We reduce these fourth-order theories to second order 
theories at the expense of introducing new scalar and tensor fields - the scalaron and the tensoron. In 
order to obtain junction conditions we apply the method of generalized Gibbons-Hawking boundary 
terms which are appended to the appropriate actions. After assuming the continuity of the scalaron 
and the tensoron on the brane, we recover junction conditions for such general brane universe models 
previously obtained by different methods. The derived junction conditions can serve studying the 
cosmological implications of the higher-order brane gravity models. 

PACS numbers: 98.80.Cq, 04.50.-h, ll.25.Mj 



I. INTRODUCTION 

In the simplest approach to the variational principle 
of any theory of gravity, such as Einstein relativity, one 
used to assume that both the variation of the metric ten- 
sor 5 g^, and the variation of the first derivative of the 
variation of the metric tensor Sg a b;c, vanish on the bound- 
ary of the integration volume [1]. However, motivated by 
quantum cosmology, Gibbons and Hawking [2] claimed 
that the latter requirement is too strong, and so they 
suggested exact cancellation of the term coming from 
the variation of the action which involves the derivatives 
with a postulated extra boundary term - nowadays called 
just the Gibbons-Hawking boundary term. For Einstein 
relativity the Gibbons-Hawking boundary term is con- 
structed of the trace of the extrinsic curvature [2] . It is 
expected that any theory of gravity may be appended 
by an appropriate boundary term. Hawking and Lutrell 
[3] found a boundary term for the fourth order gravity 
theory composed of the combination of the square of the 
Weyl tensor and an arbitrary function of the scalar curva- 
ture in four dimensions and studied its Wheeler-deWitt 
quantization. The boundary terms for gravity theory of 
an arbitrary function of the scalar curvature were also 
studied by Barrow and Madsen [4] and in a more gen- 
eral case of an arbitrary function of curvature invariants 
by Barvinsky and Solodukhin [5] . The study of the ap- 
propriate boundary terms for the Gauss-Bonnet density 
being one of the general Lovelock densities [6, 7] has 
also been done [8-11]. These densities are, however, spe- 
cial combinations of the curvature invariants which give 
the second-order field equations. The Gibbons-Hawking 
boundary term for Lovelock gravity with AdS asymp- 
totics has been found, too [12]. 

Brane gravity theories initiated by Hof ava and Wittcn 
[13] and further developed by Randall and Sundrum [14] 



'Electronic address: abalcerzOwmf.univ.szczecin.pl 
t Electronic address: mpdabfz9wmf.univ.szczecin.pl 



add an extra aspect to the problem of boundary terms, 
since the boundary of the bulk spacetime is the brane, 
which serves as the actual universe [15]. In analogy to 
a surface layer problem of clcctromagnctism, in brane 
models, one has to study appropriate junction conditions 
[16]. The problem of the gravitating surface layer junc- 
tion conditions was first solved by Israel [17], and recently 
applied to brane universes in Gauss-Codazzi formalism 
[18]. It was also solved for a Brans-Dicke braneworld [19]. 
All these considerations became the basic conditions for 
studying cosmological solutions within the framework of 
brane scenario which allows the modification of the New- 
ton's law on small scales. Brane models produce a differ- 
ent cosmological framework, since the square in energy 
density term and the dark radiation resulting from trans- 
ferring the gravitons from the bulk to the brane and vice 
versa appear in generalized cosmological equations [18]. 

Within the physically interesting context, it seems 
challenging to formulate braneworld scenario for gener- 
alized gravity theories such as Lovelock gravity and the 
fourth-order gravity. In order to achieve that, one nec- 
essarily has to formulate Israel junction conditions for 
these gravities on the brane. In fact, for the Gauss- 
Bonnet brane, these conditions were formulated by the 
application of the Gauss-Codazzi formalism in Rcf. [20]. 
This allowed many detailed studies of the Gauss-Bonnet 
brane cosmologies [21-24], and of general Lovelock cos- 
mologies [25]. However, Lovelock brane models are free 
from the problem of divergencies resulted from the ap- 
pearance of the powers of the delta function in the field 
equations [25, 26]. In fact, gravity theories which are 
based on the lagrangians being the functions of curva- 
ture invariants such as f(R) theory (see [27]; for a re- 
cent review see [28]) or / (R, R ab R ab , R abcd R abcd ) theory 
[29, 30], unavoidably lead to such divergencies and the 
formulation of the junction conditions is a non-trivial 
task. In Ref. [26] we have proposed the resolution of 
the problem for such theories by imposing more regular- 
ity onto the metric tensor at the brane position, though 
still keeping the theory to be a fourth-order. However, 
such a strong regularity of the metric at the brane may 



seem somewhat restrictive and so we have also explored 
the equivalence of these conditions to the conditions ob- 
tained for an equivalent second-order theory with an ex- 
tra scalar degree of freedom - the scalaron [31]. Wc 
have also assumed the continuity of the scalaron on the 
brane. Such a correspondence for surface layer f(R) 
universes has already been studied in Refs. [32, 33] 
and quite recently in Ref. [34] in the Gauss-Codazzi 
approach. On the other hand, in Ref. [35] the grav- 
ity theory of the linear combination f(R, R ab , R a bcd) = 
aR 2 + bR ab R ab + cR abcd R abcd (a, b, c = const.) was stud- 
ied in the Gibbons-Hawking boundary term approach. 
In these references some example cosmological solutions 
have been found. Our current task is to extend these 
considerations onto f(R 7 R a bR ab ,RabcdR abcd ) theory by 
using the Gibbons-Hawking boundary term method, also 
in the most general case, where the discontinuity of the 
new fields - the scalaron and the tensoron - at the brane 
position is allowed. 

In Section II we discuss junction conditions for the 
fourth-order f(R) gravity theory (after transforming it 
to the second order theory) by adding an appropriate 
Gibbons-Hawking boundary term, constructed of the ex- 
trinsic curvature and an extra scalar field - the scalaron. 
In Section III we derive junction conditions for a general 
f(R,R ab R ab ,R abcd R abcd ) gravity theory, first by trans- 
forming it to the second-order theory, and then by adding 
an appropriate Gibbons-Hawking boundary term, which 
is constructed of the extrinsic curvature and an extra 
rank four tensor field - the tensoron. In fact, wc bene- 
fit from the discussion of general junction conditions for 
these higher-order theories following our complementary 
approach given in Ref. [26], but we derive more general 
junction conditions - the ones which do not possess con- 
tinuity of the scalaron and the tensoron at the brane. In 
Section IV we give our conclusions. 



II. GIBBONS-HAWKING BOUNDARY TERM 

AND JUNCTION CONDITIONS FOR F(R) 

GRAVITY 

As the first example of the application of the Gibbons- 
Hawking boundary term method to derive junction con- 
ditions for brane universes, we discuss the f(R) gravity 
theory in D spacetime dimensions 
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l [ d D x^f(R) + S bulk , p , (II.l) 



where R is the Ricci scalar, \ is a constant, S bu ik, P is the 
bulk matter action, and p = 1,2. It is known [33] that 
this theory gives fourth-order field equations and that 
it can be expressed in an equivalent way by using the 
following action 



S p -- J d D xV—g{f{Q){R-Q) + f{Q)}{ll.2) 



d D - 

£>bulk,p 



where M p is the spacetime volume, Q is an extra field 
(a lagrange multiplier), K is the trace of the extrinsic 
curvature, h is the determinant of the induced metric and 
f'(Q) — df(Q)/dQ. Varying (II. 2) with respect to Q, one 
obtains an equation of motion Q = R (provided f"{Q) ^ 
0), and this is equivalent to a scalar-tensor Brans-Dickc 
gravity with Brans-Dickc parameter lobd = 0, and the 
potential V(H) = -HQ(H) + f(Q(H)) (see e.g. [31]). 
Defining H = f'(Q) (a scalaron [31]) one can rewrite 
(II. 2) in the form 



S p = I d u x^^{HR-V{H)} + S bu i h , p .(II.3) 



The variation of (II. 3) gives 



SSr, = 



M„ 



-g ab HR - HR ab (II.4) 



-g ab V(H) + H. dc (9 ab 9 cd ~ 9 {ac g b)d ) 



X rpab 

T\ bulk,p 



dV 



+ f d D - 1 xV^h2Hn d . c g^ b g^ d 6g ab 

JdM p 

+ f d D - l x^hAH, (c n d) g^ b g^ d 5g ab 



dM p 

f d D - 1 x^h{2Hg^ b g^ d n d 5g ab ), c , 

JdMp 

where n^ is a unit vector normal to dM , and it is cho- 
sen to be "outward pointing" , if spacelike, and "inward 
pointing" , if timclike (this is according to the Gauss 
integral theorem), and the induced metric is given by 
h a b = 9ab — £n a n b (e = 1 for a timclike brane, and 
e = — 1 for a spacclike brane). In our case the brane 
divides spacetime M into two parts: M\ and M<z. This 
makes convenient the usage of the symbol M p to indicate 
integral domains of (II. 4) (this is why p = 1 or p = 2). It 
should be noticed that wr. changes its direction into the 
opposite, if the calculations are done on the other side 
of the brane (that is if p changes its value from 1 to 2), 
while the vector n a remains the same. 

Let us now briefly discuss this problem in a more de- 
tailed way. At first, assume that n a is "outward point- 
ing" for p = 1. Such an assumption implies that n a is 
an "inward pointing" for p — 2. If, for p = 1, the vector 
n a is spacclike (the case of a timclike brane), then rf_ 
is an "outward pointing", and so n a = n°_, n a n a = 1. 
If, for p = 1, the vector n a is timelike (the case of a 
spacelike brane), then rf_ is "inward pointing", and so 
n a = —n°i, n a n a = 1. Similarly, if for p = 2, the vec- 
tor n a is spacelike (the case of a timelike brane), then 
n a is "outward pointing" , and so n a = — n a (because in 
the case of p = 2, n a is "inward pointing"), n a n a = — 1. 
If, for p = 2, the vector n a is timelike (the case of a 
spacelike brane), then rf_ is "inward pointing", and so 
n a = rf_, n a n a = — 1. It should also be stressed that, for 



p = 1, we have h ac rid. = e^Q, and for p = 2, we have 
h ac n d = —tKj [1]. Unifying these considerations one 
can write that 



n Ua 
h ac n d . r . 



(-l)W 
-(-1) P , 



■lyeKl 



The last term of the formula (II. 4) contains the first 
derivatives of the metric variation Sg ab . ci and in the 
Gibbons-Hawking approach [2] these derivatives should 
not necessarily be imposed to vanish. Instead, one adds 
an appropriate boundary term, which cancels the deriva- 
tives. The problem is to find its explicit form. In order 
to do so, let us first choose the foliation fulfilling the con- 
dition that n a n b . a = 0. Using this, one can derive the 
following identity 

X c . c = D e {h ea X a } + eKn a X a + eC n (n a X a ) , (II.5) 

where D e is a covariant derivative on the brane and C,-% 
is the Lie derivative in the direction of the vector field n. 
If we take X c as 



X c = 2Hg^ b g^ d n d 5g ab 



then we have 



i a X a = -(-l) p Hh ab Sg a , 



b ■ 



(II.6) 



(II.7) 



With the help of (II. 7), the formula (II. 5) can be ex- 
pressed as 



X c c = D e {h ea X a }-(-lfeKHh ab 8g ab 

- (-l) p e{C ri H}h ab Sg ab 
- (-l) p eH{C H h ab }Sg ab 

- (-l) p eHh ab {£ H 6g ab } . 



(II. 



Using the identities 



C n h ab = -2K ab , {11.9) 

h ab {C n g ab } = 2K, 

h ab {C n 5g ab } = h ab 6{Cftg ab } + eKn a n b Sg ab 

- 2D b {h b Jn a }, 

Sh ab {£ n g ab } = -2K ab Sg ab , 



hh ab 8g ab , 



and inserting them into (II. 8), after integrating out pure 
divergencies, we obtain: 



X' 



-{-l) p V^h{eh ab H, c n c 
HKn a n b + 2en b h ae D e H}Sg ab 
2(-l) p V^heKSH 
5{2(-l) p y/^heHK} , (11.10) 



The last term in (11.10) is exactly what we need to cancel 
the derivatives of the metric variation Sg ab:c . In fact, this 
is exactly the Gibbons-Hawking boundary term which 
should be appended to the action (II. 1) [3]. It reads as 



S, 



GH,p 



-2(-l)"e 



-hHKd 



D-l„ 



(11.11) 



dM„ 



Bearing in mind that 

2Hru. c g a[b 9 c]d = 

AH, {c n d) g^ b g^ d = 



-(-l) p eH(Kg ab -K ab ) , 
-2{-l) p e[g ab H, c n c - H* a n b) ] 



the variation of the full action (supplemented with the 
boundary term Sgh,p) 



S p + Sgh.p 



(11.12) 



gives 

OOs.p 



d D - 

dM p 

2n b h ea H. e + 
^S ab } Sg ab 
2(-l)*c / 



xV^h(-l) p {e [{g ab + en a n b )C H H 
HKh ab - HK ab - 2n {a W b) 



,D-X, 



-hK5H 



(11.13) 



where the bulk parts have been omitted. Full variation 
over the bulk space, separated by a brane, requires the 
variation of both of these parts separately (i.e. first for 
p = 1, and then for p = 2). This means that the full 
action is 



Op — o SiP = i + b St p = 2 + 6j 



where 



SS, 



brane 



brane •> 



d D - 1 xV^h^-S ab 6g ab 



OM„ 



(11.14) 



(11.15) 



After varying the total action (11.14), we obtain the fol- 
lowing junction conditions 

- {g ab + en a n b )[H. c n c ] - 2n {a h eb) [H. e ] (11.16) 

- [HK\h ab + [HK ab ] + 2n {a [H^] = e^S ab , 

where for any quantity A we have defined 
[A] = A + - A~ . 



(11.17) 

After some manipulations, the Eq. (11.16) can be decom- 
posed into the following set of conditions: 

[K] = , (11.18) 

S ab n a n b = , (11.19) 

S ab h ac n b = , (11.20) 

-(D - l)[H, c n c ] - D[H]K = e^S ab h ab , (11.21) 

-h ab [H, c n c ] - [H]Kh ab + [HK ab ] (11.22) 

— e ~n^° h ca h db . 



These arc the most general junction conditions for f(R) 
gravity on the brane. A physical example of a model 
which allows the discontinuity of the scalaron [H] 7^ 
on the brane is the model which possesses two different 
cosmological constants (Ai and A2) on each side of the 
brane. In such a case, the Ricci scalar is discontinuous 
and so the scalaron H = f'(R) is discontinuous, too. 

After additional assumption of the continuity of the 
scalaron at the brane position, i.e., after taking 



[H] =0 



(11.23) 



one gets less general junction conditions which were ob- 
tained in Refs. [33, 34]: 



[K] = 



[H, c n c ] = -e 



2(13-1) 



S a h ab , 



(11.24) 
(11.25) 



H[K ab ] = f ^{S cd h ca h db - 



h ab h cd scd ^ (n 26) 



III. GIBBONS-HAWKING BOUNDARY TERM 

AND JUNCTION CONDITIONS FOR F(X,Y,Z) 

GRAVITY. 

In Refs. [29, 30] a very general fourth-order gravity 
theory 



S = - I d D Xy 

XJm 



-gf(R,R ab R ab ,R abcd R abcd ) (III.l) 



has been studied. It has been shown that such a theory 
allows basic cosmological solutions such as Einstein, de- 
Sitter and Godel universes. Motivated by this result, we 
now develop similar theory, but within the framework of 
brane universes. As it has already been said, the main 
problem to start cosmological solutions is to formulate 
the appropriate Israel junction conditions for such a the- 
ory. This is what we are going to do now. 

In order to achieve the task, we apply the Gibbons- 
Hawking boundary term method to D-dimensional brane 
universes, given by the action [26] 



,5' 



1 



d D , 



X Jm 

1 ^brane < ^bulk,p 



r, f ( T? D T? ao T? r?abcd\ 



(III.2) 



jbed 



where R, R ab R , R a bcdR a arc curvature invari- 
ants, and x i s a constant. This generalized fourth- 
order theory includes the first Euler density the- 
ory: f{R,R ab R ab ,R abcd R abcd ) = xmR, as well as 
the second Euler density theory (Gauss-Bonnet term): 

ft D r> r?ab d r>abcd\ ^,., / d r>abcd a d r>ab i 

J{K,K ab K ,K abcd K ) — X K 2\ttabcdn ~iH ab H + 

R 2 ) (ki,k 2 = const.) [6], as special cases. In order to 
discuss the boundary terms for the action (III. 2), we no- 
tice that it is still a special case of a more general theory 
with the action in the form [36] 



Sg = X 



d D , 



M 



^gfigabiRabcd)- (HI. 3) 



Following the pattern of the f(R) theory given in Section 
II, we notice that (III. 2) is equivalent to the following 
second-order theory 



Si = x' 1 [ d D xy/=g{H ghij (R 

+ f(9ab,<t>cdef)} , 



•ghij Vghij , 



where 



and 



Jjghij 



_ df(g a b,(t>abcd) 



r/}2 



det 



d f{g a b, fiabed) 



VghijOpklmn 



7^0 



(III.4) 



(III.5) 



(III.6) 



The condition (III. 6) allows to get an equation of motion 



of the 



Hdi, 



field as 



'klvrin -ftghij 



R n hii in a similar manner 



as for the Q field in the previous Section II. In fact, here 
the tensor H 9hl: > generalizes the scalaron H = f'(Q) of 
(II. 2) onto a tensorial quantity. Because of that, we will 
call it a tensoron. Varying the equivalent action (III. 4), 
we get the boundary terms in the form 



X~ l { / d^xV^hAW^nadg^c (III.7) 

JdM„ 



d D - l xV^hA^ cd , d n c 5g ab } 



(>M„ 



where 



A 



abed 



\{H 



aedb 



H" 



abdc Tjcbda rjacbd 



H a 



H 



bed 1 Tjcbad 



H" 



Tjcbad\ 



(in. 



In a special case of the f(R, R ab R ab , R a bcdR abcd ) theory, 
the tensor A ahcd reads as 



A" 



bed £ t ad cb cd ba \ 

= Jx(g 9-99) 

+ f Y (2R ad g bc - R cd g ba - R ba g cd ) 

+ Af z R acbd . (III.9) 



Finally, an appropriate Gibbons-Hawking boundary term 
for the action (III. 4) is 



S, 



GH.p 



- (-1) P / 

J d Mr, 



(III. 10) 
d D ^ 1 x4^hA { - ab ^ d n c n d Lr l gab ■ 



Introducing X c = A^ ab ^ cd n d Sg ab into (II. 5), and applying 
the relations 

CfiSg ab = 5{£ftg a b} - g e b{Sn e }. a - g ea {Sn e }. b (III. 11) 



A^ cd n c n d n b = , 



(111.12) 



where the second equation (III. 12) is fulfilled for A^ ab ^ cd 
given by (III. 9), one obtains the following junction con- 
ditions for the f(R,R ab R ab ,R abcd R abcd ) brane gravity: 



[KA^ cd ]n c n d + [£ n A {ab)cd ]n c n d 

e[A (ab)cd Kcd] _ g -b [A (ef)cd Kef]ncnd 



(111.13) 



+ 2e{D s A^ cd n c n d ]h s e hfnV - 2e[A^ cd . {c ]n d) = ^S ab 



7 

n b n c [£fig ad ] - n a n c [£ n g db \ - n b n d [£ n g ac 
n a n d [£fig cb } = . 



(111.14) 



These junction conditions coincide with those obtained 
by the application of a different method in our previous 
Ref.[26] (Eqs. (5.4)-(5. 5)), provided that 



[A (Qb)cd ] = , 



(III. 15) 



i.e., after the assumption that the tensoron is continuous 
at the brane position. This is an analogous condition to 
the condition for the continuity of the scalaron on the 
brane (11.23) in the f(R) theory of brane gravity. 



IV. SUMMARY 

In this paper we derived the most general junction 
conditions for both f(R) and f(R,R ab R ab ,R abcd R abcd ) 
braneworld gravities by the application of the Gibbons- 
Hawking boundary term method. We generalized previ- 
ously obtained junction conditions for f(R) gravity for 
the case in which we did not assume the continuity of 
the scalaron field H = f'(Q) on the brane. Such a 
case appears, for example, if one takes two different val- 
ues of the cosmological constants in the bulk on each 
side of the brane. After assuming the continuity of the 
scalaron, these most general junction conditions reduce 
to the ones obtained in earlier references [33, 34]. Next, 
we derived the most general junction conditions for the 
f(R,R ab R ab ,R abcd R abcd ) braneworld gravity. Here, we 
also did not make any assumption about the continu- 
ity of the tensoron field A^ ab ' cd at the brane position. 
Again, we have shown that these junction conditions re- 
duce to those obtained by different methods in our ear- 
lier paper [26], provided one assumes the continuity of 
the tensoron field at the brane. We should stress that 
even if the scalaron and tensoron fields are discontinuous 



on the brane, and therefore their normal derivatives pro- 
duce a delta function singularity on the brane, the quan- 
tity which determines a jump of the normal derivative of 
the scalaron or the tensoron on the brane is well-defined 
(e.g. in the case of delta function having a singularity at 
some point xo, a jump of delta function xq is equal to 
zero - this comes directly form the definition of the limit 
of the function) . This obviously means that neither the 
scalaron nor the tensoron have to be continuous across 
the brane, so that the junction conditions (II.18)-(II.22) 
and (III.13)-(III.14) are physically realistic. 

We think that the method of the Gibbons-Hawking 
boundary term is more elegant than other methods of 
deriving junction conditions, provided one is able to 
suggest a correct boundary term. Up to our knowl- 
edge, so far, only the gravity theory of the linear com- 
bination of the curvature invariants f(R, R ab , Rabcd) = 
aR 2 + bR ab R ab + cR abcd R abcd (a, b, c = const.), was stud- 
ied in the Gibbons-Hawking boundary term approach 
[35]. 

We consider our result as a basic step in order to write 
down cosmological equations for f(R, R a b, Rabcd) brane 
gravity and to check, if basic csomological models are al- 
lowed in such a framework. Besides, the discontinuity of 
the scalaron and tensoron on the brane can be used to 
describe some other physical cases of the surface layers, 
such as boundary surfaces separating stars from the sur- 
rounding vacuum (singular hypcrsurfaces of higher order 
[17]). 

We also hope that the analysis of the higher-order 
brane cosmologies will give some characteristic pattern 
in statefmder (jerk, kerk/snap, lerk/crackle, merk/pop) 
diagnostic of cosmology [37], in a similar manner as it 
was given in Rcf. [38] for f(R) non- brane gravity mod- 
els. This, however, will be the matter of a separate paper 
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